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Alternative cosmologies, based on extensions ol General Relativity, predict modified thermal 
histories in the Early Universe during the pre Big Bang Nucleosynthesis (BBN) era, epoch which is 
not directly constrained by cosmological observations. When the expansion rate is enhanced with 
respect to the standard case, thermal relics typically decouple with larger relic abundances. The 
correct value of the relic abundance is therefore obtained for larger annihilation cross-sections, as 
compared to standard cosmology. A direct consequence is that indirect detection rates are enhanced. 
Extending previous analyses of ours, we derive updated astrophysical bounds on the dark matter 
annihilation cross sections and use them to constrain alternative cosmologies in the pre-BBN era. We 
also determine the characteristics of these alternative cosmologies in order to provide the correct 
value of relic abundance for a thermal relic for the (large) annihilation cross-section required to 
explain the PAMELA results on the positron fraction, therefore providing a "cosmological boost" 
solution to the dark matter interpretation of the PAMELA data. 

PACS numbers: 95.35.+d,95.36.+x,98.80.-k,04.50.-h,98.35.Gi,98.70.Sa 



I. INTRODUCTION 



Big Bang Nucleosynthesis (BBN) is the deepest avail- 
able probe of the early Universe. Its success in explain- 
ing the primordial abundances of light elements [1, 2], 
combined with Cosmic Microwave Background (CMB) 
and Large Scale Structure (LSS) studies, confirms the 
standard model of cosmology since the BBN epoch at 
MeV temperatures. At those temperatures the universe 
must have been essentially radiation-dominated. Before, 
however, a period of very enhanced expansion may have 
occured. In two earlier works [.3, 4] some of us have 
used indirect searches for Dark Matter (DM) annihila- 
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tion, namely antiprotons and 7-rays from the Galac- 
tic Centre, to derive limits on the pre-BBN expansion 
rate of the universe. The basic idea is as follows. If 
DM is composed of Weakly Interacting Massive Parti- 
cles (WIMPs) which thermalize in the Early Universe 
and then freeze-out their abundance before BBN, then 
the expansion history since freeze-out and the precise 
WMAP measurement of the dark matter relic abundance 
{nh^)^^{f^ = 0.1131 ± 0.0034 [.^] fix the annihilation 
cross-section {(Tannv) (for a given DM mass moAi)- A 
faster pre-BBN expansion requires a larger annihilation 
cross-section in order to meet the relic abundance bound, 
and in turn enhanced DM-induced astrophysical fluxes 
result. Thus it is possible to draw an upper limit on the 
Hubble rate H{T) before the BBN epoch [.i, 4]. In the 
present work we revisit this subject mainly motivated by 
the host of astrophysical data released in the last years, 
such as cosmic-ray electrons and positrons (PAMELA 
[6], ATIC [7], Fermi-LAT [n], HESS [9, 10]), antipro- 
tons (PAMELA [11]) and 7-rays (HESS [12], Fermi-LAT 
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[13, 14]). The rising behaviour of the positron fraction 
observed by PAMELA [G], in particular, has been posed 
under deep scrutiny, and in addition to astrophysical in- 
terpretations [15-17] it has been discussed the possibil- 
ity that the raise is due to DM annihilation dominantly 
occurring into leptons [18]. The DM interpretation re- 
quires large values of (annnv), which are then incompat- 
ible with a succcssfuU thermal relic. Mechanisms have 
been put forward in order to solve this problem [18, 19]. 
Furthermore, results from N-body simulations have been 
recently presented [2U, 21] allowing for a better modelling 
of dark structure and substructure in our Galaxy. 

Consequently, we are now in a position to reassess the 
constraints on pre-BBN cosmologies by using a rather 
complete scheme of observables. Besides an update 
of the previous works [3, 4], we are also interested in 
the possibility of accomodating the rising positron frac- 
tion with annihilations of a thermal DM particle whose 
properties naturally arise in non-standard cosmologies. 
Refs. [3, 4, 22-24], for instance, show that the large an- 
nihilation cross-sections required to meet the measured 
positron fraction or electron spectrum are attainable in 
the context of non-standard cosmological scenarios such 
as low reheating temperatures, scalar-tensor theories of 
Gravity, kination phases or branc world cosmology. Anal- 
yses of phenomenological consequences of modified cos- 
mologies are also discussed in Ref.s [25]. 

In this paper we pursue a generic parameterisation of 
the pre-BBN expansion rate [3, 4] and consider a multi- 
messenger, multi-wavelength scheme of DM constraints. 
We then close with an explicit simple realization of the 
cosmological enhancement mechanism in a scalar-tensor 
theory of Gravity. 



II. SURVEY OF BOUNDS ON ANNIHILATING 
DARK MATTER 

Dark Matter annihilations could be sources of elec- 
trons and positrons, (anti)protons, (anti)deuterons, pho- 
tons and (anti)neutrinos. Such particles, or their inter- 
actions in the surrounding medium, provide a plethora 
of ways to effectively constrain the intrinsic properties of 
annihilating DM (for a review see e.g. [27]). Here we will 
focus on galactic positrons, antiprotons, 7-rays and radio 
photons, and on constraints related to the optical depth 
of Cosmic Microwave Background photons. Other rele- 
vant channels, which arc not included in our analysis, are 
neutrinos from the Galactic Centre (GC) [2n] (neutrinos 
from Sun [29, 30] and Earth [30] do not directly con- 
strain the total annihilation cross section), and 7-rays 
from dwarf spheroidal galaxies [31, 32], galaxy clusters 
[33] and cxtragalactic halos [34]. For multi- messenger. 
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FIG. 1: Upper limits on the DM annihilation cross-section 
(versus the DM mass) coming from different astrophysical 
observations. Here we have considered the Via Lactea II DM 
distribution, DM annihilations into e"'"e~, the MED propa- 
gation model for cosmic rays and an electron spectral index 
7e = 3.3. The region above the thick black line is excluded 
by the convolution of all the implemented constraints. The 
shaded band labelled as "radio band" denotes the uncertainty 

on the radio constraint. The dashed line labelled as "^A, 

best-fit" denotes the values of the DM annihilation cross- 
section required to explain the PAMELA data on the positron 
fraction. The unitarity bound assuming s-wave annihilations 
[2()] is also shown. 
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FIG. 2: The same as in Fig. 1, for DM annihilations into 
W+W-. 
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FIG. 3: The same as in Fig. 1, for Aquarius DM distribution. 
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FIG. 4: Tiie same as in Fig. 2, for Aquarius DM distribution. 



multi-wavelengtli analyses see e.g. [.32, 35-38]. 

In tlie present work we adopt a model-independent ap- 
proach and consider generic annihilating DM candidates 
of masses ttidm in the interval [10 GeV, 10 TeV]. We 
study the bounds on the velocity-averaged annihilation 
cross-section (cTann^') for each annihilation channel sepa- 
rately: e+e~, /i"*"^^, T+T^, W^W~ or bb. This scheme 
basically captures the essential features of several well- 
motivated DM particles such as the lightest supersym- 
metric particle (LSP) or the lightest Kaluza-Klein par- 
ticle (LKP). We will consider throughout the paper the 
case of a DM particle which dominantly annihilates as 



FIG. 5: The same as in Fig. 1, for a cored isothermal DM 
distribution. 
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FIG. 6: The same as in Fig. 2, for a cored isothermal DM 
distribution. 



Annihilation signals are proportional to the square of 
the dark matter density, and it is therefore crucial to 
properly specify the DM distribution in the Milky Way. 
We rely on the recent results of the high-resolution N- 
body simulations Via Lactea II (VL2) [20] and Aquarius 
(Aq) [21]. In the former, the total DM profile (smooth -|- 
clumpy components) is well fitted to a Navarro-Frenk- 
White (NFW) profile, while in the latter the density 
distribution follows an Einasto profile with a = 0.17. 
Once the clumpy component has been subtracted to 
the virial mass of the simulated Milky- Way-like halos. 
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FIG. 7: Summary of the astrophysical bounds on the DM an- 
nihilation cross-section vs. the DM mass, for the Via Lactea 
II DM distribution and for different DM annihilation chan- 
nels: e"'"e~, n'^ fi~ , T^T~, W^W' and bb. Cosmic rays are 
propagated in the MED model and the electron spectral index 
is 7e = 3.3. 



FIG. 8: Summary of the astrophysical bounds on the DM an- 
nihilation cross-section vs. the DM mass, for the Aquarius 
DM distribution and for different DM annihilation channels; 
e+e~, ij.'^fi~, T^T~, W^W~ and bb. Cosmic rays are prop- 
agated in the MED model and the electron spectral index is 
7e = 3.3. 



the local density of the smooth component is Pq = 
0.41 GeV/cm^ and Pq« = 0.58 GeV/cm^ Extrapolating 
the results of the simulations down to a minimal sub- 
halo mass for WIMPs of lO"*^ Mq (e.g. [:!<)]), one gets 
^ 53(17)% of the total Milky Way mass distributed in 
virialised substructures for Via Lactea II (Aquarius) (see 
Ref. [4(1] for further details and analysis of the findings of 
the two simulations). Notice that different WIMPs and 
different early universe expansion rates may induce min- 
imal subhalo masses significantly smaller or larger than 
the fiducial value 10~^ Mq. However, this would not 
change significantly our results or conclusions. In fact, 
as shown by [40], the clump fraction becomes important 
at about 70(100) kpc for Via Lactea II (Aquarius) sim- 
ulations. This is because tidal forces deplete the pres- 
ence of subhalos in the central part of the halo. There- 
fore, the radio, 7-ray from the galactic centre (and to 
less extent, electron and positron) constraints depend 
only on regions where the subhalo component is not 
dominant. The inverse Compton scattering constraints, 
which could depend in principle on the subhalo abun- 
dances, turn out to be unimportant in our analysis, while 
the CMB constraints refer to a period where no struc- 
tures had formed thus being insensitive the problem. As 
an additional benchmark for the galactic distribution of 
DM, we consider the case of a smooth isothermal pro- 
file with no substructures. Following [41, 42], we set 
/og° = 0.4 GeV/cm^ 



The astrophysical bounds on the dark matter annihi- 
lation cross section (cTanni^)- a-s a function of the dark 
matter mass moM, for the different annihilation chan- 
nels and for the different dark matter density profiles, 
are summarized in Figs. 1-9. Due to the instrinsic com- 
plexity of these figures, we refer the reader to the inter- 
nal labellings and to the figure captions in order to pin 
down the various bounds, without excessive comments in 
the body of the article. We just remark that Figs. 1 to 6 
show the separate impact of all the astrophysical bounds, 
for specific annihilation final states and specific DM halo 
profiles. Figs. 7, 8 and 9 represent a summary of the as- 
trophysical bounds on the DM annihilation cross-section 
vs. the DM mass, for the Via Lactea II, Aquarius and 
cored isothermal halo profiles, respectively. We notice 
that the effect of the bounds may be non trivial, depend- 
ing on the DM mass, the annihilation channels and the 
halo profile: the most stringent bound arises as a combi- 
nation of different observational sources, it is tipically a 
factor of 10 (up to 100) stronger for the Via Lactea II DM 
distribution, as compared to the isothermal profile (with 
Aquarius typically staying in between). Clearly, signals 
which are more sensitive to the galactic center DM shape 
(like the radio bound or gamma-rays from the galactic 
center) induce bounds which are more affected by the 
uncertainty on the DM distribution. This is clearly seen 
in the radio band constraint in Fig.s 1, 3 and 5 for the 
e+e~ channel and Fig.s 2, 4 and 6 for the W~^W~ chan- 
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FIG. 9: Summary of the astrophysical bounds on the DM an- 
nihilation cross-section vs. the DM mass, for a cored isother- 
mal DM distribution and for different DM annihilation chan- 
nels: e"'"e~, n'^ fi" , T^T~ , W^W' and bb. Cosmic rays are 
propagated in the MED model and the electron spectral index 
is 7e = 3.3. 



ncl: moving from Via Lactca II to Aquarius, the bound 
loosens by about one order of magnitude; the bound is 
instead vifeakened by about 4 orders of magnitude for an 
isothermal DM distribution. Other signals, like e.g. an- 
tiprotons, are not very sensitive to the DM halo profile 
[4-5], as can be seen by confronting again Fig.s 2, 4 and 
6. Instead, signals like antiprotons or positrons are more 
sensitive to the diffusive region of the galactic environ- 
ment [44, 45], and therefore will produce bounds which 
are limited by the uncertainties on the cosmic-rays prop- 
agation properties: the effect of this will be shown in 
Section IV. 



A. Antimatter 

Unlike neutral particles, positrons and antiprotons pro- 
duced in the Milky Way undergo different processes 
that change their direction and energy while crossing 
the galactic medium. The galactic magnetic fields, 
for instance, are responsible for deflection and, due to 
their (poorly known) inhomogeneities, the evolution of 
a positron or an antiproton can be treated as a random 
walk with a certain diffusion coefficient. For positrons, 
important phenomena to take into account are energy 
losses through inverse Compton scattering off the CMB 
and starlight and synchrotron emission, which proceed 
at a space-independent rate b{E^+) ~ i?^^./(GeV • te) 



with Te — 10"'^^ s [4.3, 44] and result in a diffusion length 
of a few kpc. Antiprotons, on the other hand, can 
travel larger distances without losing much energy by 
synchrotron or inverse Compton since nip 3> m^. In- 
stead, they may be swept away by galactic winds, as- 
sumed to be constant and perpendicular to the disk. Fur- 
thermore, annihilations pp mainly in the disk are respon- 
sible for the disappearance of primary antiprotons. To 
compute propagated fluxes wc apply the formalism out- 
lined in [45, 46] (and references therein) and use the MIN, 
MED and MAX propagation models. In this framework, 
the total positron or antiproton flux at Earth for a spe- 
cific dark matter candidate and a certain propagation 
model is the sum of the smooth and clumpy DM compo- 
nents: 



b^^{E) = {l-fQf4>UE) + {4>ci){E) 



(1) 



where /© = 3.0(0.23)% is the local clump fraction for 
VL2 (Aq) and E plays the role of kinetic energy in the 
antiproton case. Recall that both (ps-m and (0ci) are pro- 
portional to ((Tannw). We notice that the term has 
been modelled following [46] and folds self-consistently 
the population of subhalos and the Green functions for 
the propagated positrons/antiprotons. 



1. Electrons and positrons 

The PAMELA satellite has measured the positron frac- 
tion (/)e+/((/)g+ 4- (/)g-) in the energy range 1.5—100 GeV 
[()]. The data show a steep rise above ^ 7 GeV. Here 
we restrict ourselves to i?g+ > 10 GeV and therefore 
disregard solar modulation. In order to translate the 
PAMELA positron data into an upper limit on the DM- 
induced positron flux (t)^+^ ^ we assume (j) a power-law 
electron flux 0e- oc E~J'' , 7e = {3.3,3.4,3.5} (in rough 
agreement with [47]) normalised to the AMS-01 mea- 
surement at 10 GeV [4S]; and (ii) a secondary positron 
flux (f)'^^'^ produced by the spallation of cosmic rays in the 
Galaxy [43] . In this scheme, 0e+ = 4'f+^ + 4'^'+ ■ We now 
require, for fixed rriDM, that no single energy bin of the 
PAMELA positron data above 10 GeV is exceeded by 
more than 3a. This produces an overall upper limit on 



(fann^')- A Complementary bound, particularly for multi- 
TeV leptophilic DM candidates, is motivated by the re- 
cent measurements of the electron plus positron flux by 
Fcrmi-LAT [s] and HESS [9, 10]. As with PAMELA 
positron fraction, we disregard solar modulation - Fermi- 
LAT and HESS cover energies ranging from ~ 25 to ~ 
5000 GeV - and, conservatively, draw the 2a upper limit 
on ((Tanni') by Only considering electrons and positrons 
produced by DM annihilations. 
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On the other hand, we are also interested in study- 
ing the feasibihty of explaining the PAMELA positron 
fraction with DM annihilations. So, for fixed nioM, we 
fit the data above 10 GeV leaving (aannv) as a free 
parameter and demand that < 20 (7 data points 
are available). Whenever the positron best-fit values 
of (o'ann''^) are not in conflict with bounds coming from 
other astrophysical signals, we consider that value as a 
positive solution of the PAMELA data in terms of dark 
matter annihilation. More detailed fitting procedures 
have been presented in other references e.g. [49]. It is 
well-known that it is possible to obtain good fits to the 
PAMELA positron data with DM particles which annihi- 
late preferably into leptons (hadrons) with masses above 
~ 100 GeV (a few TeV). Of course, the required annihi- 
lation cross-sections are larger than the standard thermal 
value, 2.1 • 10"^^ cni^s"^. Wc clearly recover these results 
in our analysis. Differences with respect to Refs. such as 
[49] are due to a different local dark matter density (we 
are using values greater than the "usual" 0.3 GeV/cm'^), 
the inclusion of dark matter substructure (according to 
Via Lactea II and Aquarius simulations), and the use of 
a slighlty smaller electron fiux (we normalise to AMS-01 
at 10 GeV and consider a power-law index 3.3 as a refer- 
ence value). All such factors play to lower our "PAMELA 
best-fit" cross-sections with respect to other Refs. 

Note as well that a joint explanation of the PAMELA 
positron fraction and the Fermi-LAT/HESS electron 
plus positron fiux in terms of DM annihilations favours 
heavy rather than light DM particles. However, in this 
work we do not pursue a global fit to PAMELA/Fermi- 
LAT/HESS but instead use the electron plus positron 
fiux as an astrophysical constraint only. 



2. Antiprotons 

The comparison of the PAMELA antiproton ratio 
4'p/4'p [11] with theoretical estimates of secondary an- 
tiprotons reveals little space for p deriving from DM an- 
nihilations (or any other primary source) [30] . These con- 
siderations disfavour light DM particles decaying promi- 
nently into hadrons. 

To derive the antiproton bound we use the interstellar 
proton flux <f)p and the interstellar secondary antiproton 
flux (p^p"^ as given in Ref. [>!()], and apply a solar mod- 
ulation in the force field approximation with 4>f = 500 
MV. The 2(7 upper bound on (cTann^) from antiproton 
searches is derived by using the whole energy spectrum 
in the range where PAMELA antiproton measurements 
are available [11]. 



B. 7-rays 

In some specific models, DM particles can annihilate 
directly into photon(s) producing clear spectral 7-ray 
lines (see e.g. [olJ]): the branching ratios for such pro- 
cesses are usually rather low. We do not consider in this 
paper such annihilation channels. Instead, we consider 
7-rays produced by decays of (or radiation from) final 
state particles of the annihilation process. These pro- 
cesses lead to a continuous spectrum up to energies close 
to the mass of the DM particle (e.g. [-"jl, 52]): 

where dN^/dEj is the differential spectrum per annihi- 
lation for a given annihilation channel [37, 53] (we re- 
mind that we assume all the annihilation proceeds via a 
single channel). We compute the 7-ray fiux in a 10~^ 
sr solid angle towards the Galactic Centre (GC). We 
compare our predictions with the HESS measurement 
of the GC source in 2003 and 2004, $(> 160 GeV) = 
1.89±0.38 X 10"" 7 cm"2g-i ^i2j^ and derive the bound 
on the annihilation cross-section requiring that our model 
does not exceed the measurement by more than 2a: this 
is the bound from the galactic center, labelled as "7 from 
GC" in Figs. 1-6. 

Another way for DM annihilations to give rise to 7- 
ray fluxes is through inverse Compton scattering (ICS) 
on electrons and positrons produced by DM annihilation. 
In fact, low-energy photons, such as those in the CMB, 
starlight and infrared radiation, may be up-scattered by 
high-energy electrons and positrons. This channel has 
gained particular relevance ever since the PAMELA col- 
laboration published their exciting results: if the positron 
data arc due to DM annihilations, then there should exist 
a large population of electrons and positrons in the Milky 
Way able to up-scatter low-energy photons. In Ref. [55] 
the authors computed the ICS 7-ray spectrum towards 
regions far from the GC. Here we apply the same proce- 
dure. Notice that these constraints are especially robust 
against the DM density profile since the GC region is 
excluded from the field of view. 



C. Radio photons 

The interpretation of the rising positron fraction ob- 
served by PAMELA as a signal of DM annihilations 
would result in a large amount of highly energetic elec- 
trons and positrons permeating our Galaxy. This should 
hold true particularly towards the GC where the DM 
density is expected to be the highest. Such relativistic 
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electrons and positrons propagating in the galactic mag- 
netic field emit synchrotron radiation in the radio fre- 
quency band. Let us focus on a region towards the GC, 
small enough so that diffusion does not play an impor- 
tant role and where the galactic magnetic field is strong 
enough to neglect electron (and positron) energy losses 
other than synchrotron emission. Assuming further that 
advection is negligible and following Refs. [32, 37], the 
total synchrotron power emitted by the distribution of 
DM-induced electrons and positrons is: 



syn 



2 

, (3) 

where Ep{y.,v) = Airralv / {?, ■ 0.29eB(x)). We com- 
pute solely the contribution given by the smooth distribu- 
tion of DM. Following Ref. [')(>], we consider a cone with 
half-aperture 4" pointed toward the GC and ly ~ 0.408 
GHz for which an upper bound of 0.05 Jy has been de- 
rived from radio observations [.^7]. In this way we can 
determine the constraint on the plane moM — (fTannw) 
given by radio observations of the GG. In order to take 
into account the different approaches of Refs. [3")] and 
[32], we weaken the bound obtained with Eq. (3) by a 
factor 7 [77] and consider such rescaled result as an ef- 
fective constraint. The radio bound is shown as a shaded 
band in Figs. 1-6. 



D. Optical depth of CMB photons 

The optical depth of GMB photons depends on how 
and when the re-ionization of the universe occured. If 
DM is annihilating, then a considerable amount of high- 
energy electrons and positrons may be created after re- 
combination giving rise, by ICS on CMB photons, to a 
population of (low-energy) 7-rays. These can easily ion- 
ize the gas releasing electrons and hence reducing the op- 
tical depth of GMB photons [-jS]. Comparing such effect 
with the measured optical depth by WMAP, the authors 
of Ref. [59] have derived upper limits on the annihilation 
cross-section of DM particles, that we also reproduce in 
Figs. 1-6. Since there is clearly no dependence on the 
DM density profile, these constraints are very robust and 
difficult to avoid. This constraint has also been carefully 
analyzed in Refs. [GO] and [fil]. 

Other effects of conspicuous DM annihilations in 
the early Universe are the heating of the intergalactic 
medium [59] and the distortion of CMB anisotropics and 
polarisation [02]. 



III. MODIFIED COSMOLOGIES 

Cosmological models arising in modification of General 
Relativity (GR) very often predict a cosmological history 
with an expansion rate H{T) larger, at early times, than 
the Hubble expansion rate Hq^{T) of standard cosmol- 
ogy. Generically, we can encode this enhancement into a 
temperature-dependent function A{T) as [3, 4, 22]: 



H{T) = A{T)Hgr{T) 



(4) 



with A{T) > 1 at large temperatures and with A{T) 1 
before BBN sets up, in order not to spoil the successful 
predictions of BBN on the abundance of primordial light 
elements. 

With an increased expansion rate, thermal relics 
freeze-out their abundance earlier than in standard 
cosmology: this implies that a thermal (cold) relic 
matches the correct relic abundance for annihilation 
cross-sections ((Tann^^) which are larger than in standard 
cosmology. A consequence of this is that dark matter par- 
ticles possess today, in the galactic environment, larger 
annihilation cross-sections and thence enhanced indirect 
detection signals, as compared to those obtained for a 
thermal decoupling in GR. This implies that indirect 
searches for dark matter, like those discussed in the pre- 
vious Section, may have a potential of constraining pre- 
BBN cosmological histories, under the assumption that 
dark matter is a thermal relic. 

We have investigated the consequences of these phe- 
nomena in Refs. [3, 4], where we used the cosmic-ray 
antiproton and gamma-ray data to derive bounds on the 
admissible enhancement of the expansion rate in the pre- 
BBN phase. In the current paper we extend these anal- 
yses to comprehend all the astrophysical observables di- 
cussed in the previous Section. 

In the case of positrons, the recent measurements from 
PAMELA on the positron fraction have shown a steady 
"anomalous" rise at energies above 10 GeV, up to about 
100 GeV (the current largest probed energy). The in- 
terpretation of this so-called PAMELA excess is cur- 
rently under deep investigation. Astrophysical origins 
of this behavior have been shown to be able to ex- 
plain the PAMELA data: local sources, like pulsars [15], 
or positron production mechanisms occurring inside the 
sources of cosmic-rays [16, 17] are suitable to reproduce 
quite well the PAMELA result. An alternative solution is 
offered by dark matter annihilation: in this case, it has 
been clearly shown in many independent analyses that 
in order to explain the PAMELA data, the dark mat- 
ter candidate needs to meet a number of requirements: 
first of all, the size of the annihilation rate has to be 
orders of magnitude larger than the one obtained for a 
thermal relic in standard cosmology i.e. referring to an 
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annihilation cross-section of (era 



2.1 • 10" 



s ^, which is the one required to obtain the correct rehc 
abundance for cold DM (17/i2)WMAP ^ o.U in standard 
cosmology; secondly it has to dominantly annihilate into 
Icptons, unless it is quite heavy (above the TeV scale), in 
this latter case it is allowed to decay also into hadronic 
channels, where antiprotons are produced. Since alter- 
native cosmologies with A{T) > 1 imply that the cor- 
rect relic abundance of a relic particle is obtained with 
larger annihilaton cross-sections, they offer a framework 
to explain the PAMELA data without requiring specific 
mechanisms to boost the annihilation rate like, for in- 
stance, Sommerfeld enhancements [18, 19] or (unlikely) 
large astrophysical boosts [4(>]. In a sense, alternative 
cosmologies offer a "cosmological boost" to a thermal 
relic. In the first part of the next Section we will discuss 
under what conditions PAMELA data are explained by 
means of this cosmological boost. 

In order to be as general as possible, we will perform 
our analysis by parameterizing the temperature-behavior 
of the enhancement function A{T) as: 



tanh 



(5) 



for temperatures T > and A{T) ^ 1 for T <T,c. This 
form has been adopted in our previous analyses [■:(, 22]: 
it is a suitable parameterization to describe a cosmology 
where H — )- i?GRi at some "re-entering" temperature T^e. 
We must require T^c > 1 MeV to make sure not to be in 
conflict with the predictions of BBN. For definiteness we 
will fix Tro = 1 MeV in our analysis, except when ex- 
plicitely mentioned otherwise. Notice that a sensitivity 
on this parameter is expected, as discussed in Ref. [22] 
and as it will be shown in the next Section. The maxi- 
mal enchancement on the relic abundance is obtained for 
the lowest possible value of T^c = Tbbn- By fixing this 
parameter at 1 MeV, we derive the lowest bounds on the 
enhacement parameter 77: larger 77 are expected for larger 

For T > Tie we have: 



^(r) = i + r/( I 



T 



(6) 



where the last implication is valid for rj':$> 1. Thus, Tf is 
the normalization temperature at which A{Ti) = 1 -(- ry, 
which, again for large values of 77, goes to A{Ti) — > 77. 
For definiteness, and to conform to our previous analy- 
ses, we take Tf to be the temperature at which the WIMP 
DM candidate freezes out in standard cosmology (Tf is 
therefore DM mass-dependent). Therefore 1 -I- 77 repre- 
sents the enhancement of the Hubble rate at the time 
of the WIMP freeze-out. The freeze-out temperature is 



determined with the standard procedure, which can be 
found for instance in Refs. [.3, 22]. 

We will organize our discussion in terms of bounds on 
77 for different cosmological models, characterized by the 
temperature-evolutionary parameter v: v = 2 refers to 
the Hubble rate evolution in a Randall-Sundrum type II 
brane cosmology scenario of Ref. [ ] ; = 1 is the typical 
kination evolution, discussed e.g. in Ref. [G4]; 1/ = — 1 is 
representative of the behavior found in scalar-tensor cos- 
mologies in Ref. [22], to which we refer for additional dis- 
cussions on thermal relics in cosmologies with enhanced 
Hubble rate. The trivial case v = Q refers to an overall 
boost of the Hubble rate, like in the case of a large num- 
ber of additional rclativistic degrees of freedom in the 
thermal plasma. 

In Section V we will instead present a specific cosmo- 
logical model arising in scalar-tensor theories of gravity, 
where we will explicitely show the ability of such mod- 
els to provide the right amount of cosmological boost to 
explain the PAMELA results. 

For more details on the modified cosmological scenar- 
ios, the calculation of the relic abundance in these mod- 
els, including some analytical results and discussion, we 
refer to Refs. [3, 4, 22]. 



IV. COSMOLOGICAL BOOST: PAMELA AND 
BOUNDS ON MODIFIED COSMOLOGIES 

We start our analysis by studying under what condi- 
tions the modified cosmological scenarios can explain the 
PAMELA data without violating any of the constraints 
described in Section II. The study is performed as fol- 
lows: we numerically solve the Boltzmann equation for 
the evolution of the number density of a thermal relic in 
a cosmology where the Hubble rate is given by Eqs. (4,5) 
and we determine the values of the enhancement parame- 
ter 77 which are required, for any given cosmology labelled 
by the parameter in order to reproduce the correct 
value Q.^h'^ = {^h'^)c^M^ = O-H for tl^e rehc abundance 
of the dark matter particle. The annihilation cross sec- 
tions arc fixed to explain the PAMELA data (within the 
implemented astrophysical bounds), and have been de- 
rived in Section II. The analysis is performed separately 
for each of the possibile DM annihilation channels. An- 
alytic considerations may be found in Ref. [■]]. 

First of all, we show the effect of the modified ex- 
pansion rate on the dark matter decoupling: this will 
be useful to understand the following analysis. Fig. 10 
shows the temperature evolution of the comoving abun- 
dance Y = n/s, where n is the number density and 
s the total entropy density. The almost-vertical line 
refers to the Boltzmann-suppressed equilibrium abun- 
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FIG. 10: Temperature evolution of the comoving abundance 
Y. The almost-vertical solid line shows the equilibrium abun- 
dance, while the lowest solid line which asymptotizes to the 
relic abundance value refers to the solution of the Boltzmann 
equation in standard cosmology (GR). The other solid lines 
refer to the solutions of the Boltzmann equation for differ- 
ent cosmologies with T^a = 1 MeV: from top to bottom 
v = -1,0,1,2. Dot-dashed lines refer to T^-c = 100 MeV 
and dashed lines to Trc = 1 GeV. All the lines are normal- 
ized to the asymptotic value of the abundance in standard 
cosmology. The mass of the dark matter particle and annihi- 
lation cross section have been fixed at: — 100 GeV and 
(a^nnv) = 2.1 ■ 10-2" cm^ s-\ 

dance. The other lines refer to the solution of the Boltz- 
mann equation for various cosmologies, and explicitely 
show the decoupling and asymptotization of Y, from 
which the relic abundance is determined as usual as: 
Q,~^h? = rriy^soYo/ pc, where sq and Yq = Y{T = Tq) arc 
the current values of the entropy density and of y, and 
Pc denotes the critical density of the Universe. All the 
curves are normalized to the asymptotic value obtained 
in standard cosmology (GR, for General Relativity), in 
order to explicitely show the amount of the enhancement 
obtained in the different cosmological models. The mass 
of the dark matter particle and annihilation cross section 
have been fixed, in this example, at: rriy, = 100 GeV 
and (aannv) = 2.1 • 10~^^ cm^ s~^. The solid lines refer 
to cosmologies labelled by the value of the v parameter, 
enhancement factor 77 = 100 and re-enetering tempera- 
ture Tre = 1 MeV: from top to bottom, 1/ = —1,0, 1,2. 



FIG. 11: Values of the i] parameter required to explain the 
PAMELA data together with Q^h^ = (n/i^)^^^, in the 
case of DM annihilation into e^e~ and for different cosmolo- 
gies labelled by the values of u. The DM annihilation cross 
section is required to explain the PAMELA data without vi- 
olating the astrophysical bounds (see Figs. 1, 3 and 5 for this 
case) and the values of ry are determined, for each DM mass, 
in order to have the correct DM relic abundance in the modi- 
fied cosmology, which therefore produces the required cosmo- 
logical boost. The solid (red) lines refer to the Via Lactea 
II DM distribution, the dot-dashed (black) lines refer to a 
cored isothermal sphere and the dotted (blue) lines refer to 
the Aquarius DM distribution. Propagation parameters are 
set at the MED case. 



We notice that larger values of v endow smaller enhance- 
ments of the relic abundance, for the same value of 77: 
this is due to the fact that a larger v implies an Hub- 
ble rate which evolves faster in temperature, and in this 
case H{T) paces closer to the fast-decreasing annihila- 
tion rate, inducing a longer phase of post-freeze-out an- 
nihilation, which in turn reduces the abundance. This 
fact implies that larger enhancement factors are required 
for large values of v in order to boost the relic abundace 
to the WMAP value, when starting with a very large an- 
nihilation cross section, like the case for explaining the 
PAMELA data. We will find this behaviour in the fol- 
lowing analysis. Fig. 10 shows also the effect of anticipat- 
ing the temperature Tj-c at which the modified expansion 
rate recovers the GR behaviour: dot-dashed lines refer 
to Tro = 100 MeV, while dashed line to T,^ = 1 GeV. We 
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FIG. 12: The same as in Fig. 11, but showing the uncertainty 
band coming from cosmic-ray propagation parameters: for 
each case, the upper hne refers to the MIN model, the lower 
line to the MAX model. Only cosmologies with v = and 
V = 2 are shown. The DM distribution is a cored isothermal 
sphere. 



notice, in these cases, that early re-entering phases are 
accompanied by a reduction in the rehc abudance after 
Tic. This phenomenon was explained in Ref. [22], where 
it was named "re-annihilation" : the (sharp) drop in the 
enhanced Hubble rate at Tj-c may result in a new efficient 
phase of annihilation of the relic particles, because their 
annihilation rate returns to be larger than the expan- 
sion rate for some time. This effect is clearly manifest 
in Fig. 10: in this case, larger enhancement parameters 
77 would be required in order to boost enough the relic 
abundance to the WMAP value, when (cTann^^) is large. 
We remind again that we will always use Tro = 1 MeV, 
unless explicitely stated. Finally, we wish to comment 
that also the speed at which the enhanced Hubble rate 
recovers GR may modify the final relic abundance. We 
modelled this phase with the tanh behaviour of Eq. (5): a 
change in this behaviour may (relatively slightly) change 
the results shown in Fig. 10. For definiteness, we will use 
the form of Eq. (5) throughout the Paper. 




m^(GeV) 



FIG. 13: The same as in Fig. 11, for DM annihilation into 
fJ^^ l-i^ (except for the Aquarius case, not shown here). The 
dotted lines refer instead to the case of Trc ~ 100 MeV, for 
the cored isothermal dark matter distribution. 

A. Boosts for PAMELA 

Let us now turn to the discussion of cosmological mod- 
els featuring a thermal relic able to reproduce the right 
amount of positrons to explain the PAMELA data with- 
out being at odds with other astrophysical data, and pre- 
senting the correct relic abundance. This is an alterna- 
tive solution to the PAMELA "anomaly" in terms of dark 
matter annihilation: the compatibility between the large 
annihilation cross-sections required by the PAMELA 
data and the WMAP value of the relic abundance are 
obtained by means of modified cosmologies. 

Fig. 11 refers to the case of annihilation into an e~^e~ 
final state, and the annihilation cross-sections used to 
obtain Fig. 11 are those shown in Figs. 1, 3 and 5 
which refer to the " _ best-fit" curve restricted to 

e+ -\-e 

the range where it is not excluded by any of the con- 
sidered constraints. The values of r] required to explain 
the PAMELA positron fraction (and compatible with the 
other bounds) change significantly with the cosmological 
scenario: they are confined in the range between 0.1 and 
100 for 1/ = — 1, while for 1^ — 2 they are significantly 
larger, being in an interval from 0(1) to about 10"^, de- 
pending on the dark matter mass and on the dark matter 
halo profile. The values of rj increase with because 
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FIG. 14: The same as in Fig. 11, for DM annihilation into 
r+r^ and into W^W' . 

the PAMELA data require larger annihilation cross sec- 
tions for larger masses, as seen in Fig. 1: in this case, in 
order to match the WMAP value for the relic abundance, 
larger cross sections require earlier decoupling of the dark 
matter particle, and this in turn requires a faster expan- 
sion of the Universe, hence larger 77. In addition, from 
Fig. 11 we notice that the enhancement factors are signif- 
icantly larger for larger values of i^, in accordance to the 
previous discussion in relation with Fig. 10. Fig. 11 also 
shows that the uncertainty arising from different choices 
of the dark matter distribution in the halo may be rele- 
vant: the solid lines refer to the case of the Via Lactea II 
dark matter distribution, the dot-dashed lines to a cored 
isothermal sphere and the dotted lines to the Aquarius 
simulation. The results shown in Fig. 11 have been ob- 
tained for the MED set of astrophysical parameters gov- 
erning cosmic-ray diffusion. 

The effect induced by the uncertainties on galactic dif- 
fusion is shown in Fig. 12: the yellow band encompasses 
the variability on the required values of 77, when the 
propagation parameters are changed from the MIN (up- 
per line of each set of curves) to the MAX (lower lines) 
values. The range of masses differs, when changing the 
propagation parameters, because the compatibility of the 
explanation of PAMELA data with the other bounds on 
indirect searches pins down different mass intervals, as 
discussed in Section II. 



FIG. 15: Upper bound on the 77 parameter for different cos- 
mological models and for the case of DM annihilation into 
e^e~. The bounds arise from the astrophysical constraints on 
the DM annihilation cross-section and from the requirement 
that the DM relic abundance matches the WMAP value for 
CDM. The solid (blue) lines refer to a cored DM distribution, 
the dashed (black) lines to the Via Lactea II DM distribu- 
tion and the dotted (red) lines to the Aquarius DM profile. 
Propagation parameters are set at the MED case. The set 
of solid, dashed and dotted lines refers to cosmologies with 
= — 1, 0, 1, 2 going from the lower to the upper curves. For 
each cosmology, the excluded values of 77 are those above the 
corresponding line. The red dot-dashed lines refer to a cored 
isothermal DM distribution and the MIN set of propagation 
parameters. 



Figs. 13 and 14 report the enhancement factors re- 
quired in the case of annihilation into fi~^fi~ ^ t~^t~ and 
W^W~ final states. In these cases, the propagation pa- 
rameters are set at the MED configuration. The solid, 
dot-dashed and dotted lines again refer to Via Lactea 
II, cored isothermal sphere and Aquarius dark matter 
distributions, respectively. The range of masses able to 
explain PAMELA data and compatible with astrophysi- 
cal observations are obtained by the analysis of Section 
II: we notice the significant difference for the T+r~ case, 
where the isothermal distribution is compatible with the 
PAMELA data for a mass interval much larger than for 
the Via Lactea II or Aquarius cases. For the W^W~ 
final state. Via Lactea II provides only marginal com- 
patibility with the PAMELA data, and is therefore not 
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FIG. 16: The same as in Fig. 15, for DM annihilation into 

present in Fig. 14. Fig. 13, for the fi^ fJ-~ channel, shows 
also the effect of anticipating the recovering of the GR 
cosmic evolution: the dotted lines refer here to Tre = 100 
MeV, for the isothermal dark matter distribution. As 
discussed above, in this case the predicted relic abun- 
dance is lower than for the case of lower values of Tic- 
as a consequence, a larger enhancement rj is necessary to 
match the correct value of the relic abudancc. This effect 
similarly affects all the other results shown in Figs. 11, 
12 and 14 and it will affect also similarly the results on 
the bounds, to which we now turn. 

B. Astrophysical bounds on modified cosmologies 

The astrophysical bounds on the dark matter annihi- 
lation cross-section discussed in Section II may be al- 
ternatively used to set constraints on the cosmological 
histories, as it was done in Refs. [3, 4], where we used 
antiproton and gamma-ray data and a dark matter par- 
ticle annihilating dominantly into a quark-antiquark fi- 
nal state (namely bb). We now extend that analyses by 
considering the whole host of experimental data of Sec- 
tion II and by including the whole set of annihilation final 
states of a cold dark matter particle. 

The results are shown in Fig. 15 for the e+e~ anni- 
hilation channel, in Fig. 16 for the fJ-^fJ.'' channel, in 
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FIG. 17: The same as in Fig. 15, for DM annihilation into 



Fig. 17 for the r+r^ channel, in Fig. 18 for W^W~ 
and in Fig. 19 for the annihilation into a bb pair. In 
every figure, the solid lines refer to a cored dark mat- 
ter distribution, the dashed lines to the Via Lactca II 
numerical simulation and the dotted lines to the Aquar- 
ius DM distribution. Propagation parameters are set at 
the MED configuration. The bound for each cosmology 
(y = —1,0,1,2, from bottom to top) is the area above 
the corresponding line. We notice that, depending on the 
dark matter mass and on the annihilation channel, the 
bounds may be quite restrictive. This is an interesting 
result, since it imposes strong bounds on the cosmolog- 
ical histories of the Universe at the time of dark matter 
freeze-out (from T - 400 MeV to T - 400 GcV for the 
mass range considered here), under the hypothesis that 
dark matter is a thermal relic. The bounds are typi- 
cally stronger for lighter dark matter, since for lighter 
dark matter the astrophysical bounds are stronger. The 
dependence on the dark matter profile may be large and 
the size of the difference depends also on the annihilation 
final state of the DM particle. The same (a dependence 
on the final annihilation state) also occurs in determin- 
ing whether a bound is stronger for Via Lactca II or for 
Aquarius (the cored isothermal halo being always less 
constraining): this occurs because the absolute bound 
on the DM annihilation cross section has origin from dif- 
ferent signals for different masses and the impact of the 
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FIG. 18: The same as in Fig. 15, for DM annihilation into 
W+W-. 

DM halo profile affects differently the various signals, as 
can be seen for the analysis of Section II. 

The effect of the propagation parameters of charged 
cosmic-rays is also shown in Figs. 15 to 19, where the 
dot-dashed lines refer to the MIN configuration. We no- 
tice that for the MIN case the bounds are much looser 
in Figs. 18 and 19: this is easily understandable from 
the fact that these are hadronic annihilation channels 
and, especially for light dark matter, the bound comes 
from antiproton searches. The propagation parameters 
induce a large uncertainty on the antiproton flux [4-5], 
and the MIN configuration predicts almost an order of 
magnitude smaller flux than the MED case: this implies 
that compatibility of the antiproton flux with the data 
allows larger annihilation cross-sections, and therefore 
looser bounds on are obtained. 
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FIG. 19: The same as in Fig. 1.5, for DM annihilation into bh. 
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FIG. 20: Evolution of the (squared) enhancememt factor of 
the Hubble rate in our modified cosmology, as a function of 
the temperature of the Universe. The solid line corresponds 
to the scalar field parameters <^in = 1.9, v'l'n = 0.45 and 6 = 8 
(Model 1); the dashed line to tp^ri ~ 1-5, (/sj^ — 0.4 and 6 = 8 
(Model 2); the dot-dashed line stands for (pin = 1.5, ip'i^-^ = 0.4 
and 6 = 4 (Model 3). 



V. A SIMPLE MODEL 

We mentioned in the previous Section that alternative 
cosmologies, based on extensions of General Relativity, 
can provide the large boost factors required to explain in 
terms of a thermal relic the positron fraction measured by 
PAMELA. We now focus on a simple model which allows 
an explicit calculation of such cosmological boosts. 



We consider the action: 

S = Sg + Sm , (7) 

where: 

(8) 
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FIG. 21: Evolution of the scalar field as a function of the 
temperature of the Universe. The solid line corresponds to 
ipiu = 1.9, v^in = 0.45 and 6 = 8 (Model 1); the dashed line 
to (pin = 1.5, <p(n = 0.4 and 6 = 8 (Model 2); the dot-dashed 
line to ((9in = 1.5, = 0.4 and 6 = 4 (Model 3). 
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FIG. 22: Factor (1 — Siu) as a function of the temperature of 
the Universe. Notice the fiuctuations associated to the differ- 
ent particles becoming non-relativistic in the thermal bath. 
In this plot we considered a "SUSY-like" spectrum. The 
dashed, solid and dot-dashed lines refer to DM candidates 
with a mass of 50 GeV, 200 GeV and 1 TeV, respectively. 

describes the gravitational interaction, which is now me- 
diated by the metric g^j^y and a scalar field ip. Sm repre- 
sents the action for the (visible and dark) matter sector. 
We assume an universal metric coupling between gravity 
and matter, that is: 

5™ = 5™[*™,^2((p).g^.] , (9) 

with 5'j„ indicating a generic field of the matter sector 
coupled to the metric A^{(p)g^u- In the present discussion 
the potential V{ip) plays no role and we therefore set it 
to zero. Thus, the function A{ip) univocally specifies the 
model. 

If A{ip) is constant, then the action Sg + Sm is just that 
of General Relativity, plus a minimally coupled scalar 
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FIG. 23: Relic abundance as a function of the DM mass, cal- 
culated for the three cosmological models of Fig. 20 and for a 
DM annihilation cross section able to explain the PAMELA 
data within the astrophysical constraints for annihilation into 
e^e~, as shown in Fig. 1. The dashed line shows the cor- 
respoding relic abundance in standard cosmology. The DM 
distribution is from Via Lactea II and the cosmic-ray prop- 
agation parameters are set at the MED case. The horizon- 
tal thin band shows the WMAP range for CDM abundance, 
inh^)cDM'' = 0.1131 ± 0.0034 [r>]. 

field. The scalar contribution to the gravitational inter- 
action is therefore measured by the quantity: 

dlogA{ip) 

If a{f) 0, then Eq. (7) defines a Scalar- Tensor theory 
of gravity formulated in the Einstein frame [Go]. 

We now consider an homogeneous cosmological space- 
time: 

ds^ ^df ~a^{t) df , (11) 

where the matter energy-momentum tensor, T^j^ = 
2{—g)~^^'^SSm/Sg^'^ admits the perfect-fluid representa- 
tion: 

T^iy ^ [p + p) u^u^ -pg^u , (12) 

with g^n u^u'^ = 1. 

The Friedmann-Robertson- Walker (FRW) equations 
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FIG. 24: The same as in Fig. 23, for DM annihilation into 



FIG. 25: The same as in Fig. 23, for DM annihilation into 



then take the form: 



[p + 3p + 2Ml^ip^] 



a 

(p + 6-(p 

a 



A/2 



1 



[a{p - 3p)] 



where p ~ w p and the Bianchi identity reads: 
d{pa^) + p da^ = {p - 3 p) a^d\ogA{ip). 



(13) 



(14) 



In the previous Sections we calculated the particle 
physics rates and cross sections in the Jordan frame, 
which is defined by the Weyl rescaling gj^ = A^{(p)g^i,„ 
(see Eq. (9)) [66]. For consistency, we thus study the 
evolution of the Universe in the Jordan frame. The new 
degree of freedom ip modifies the General Relativity Hub- 
ble expansion as follows: 



1 - (^')V6 



(15) 



where -ffcR is the General Relativity Hubble rate and a 
prime denotes a derivative with respect to the logarithm 
of the scale factor. 

We now focus on the time (or temperature) evolution 
of Eq. (15) in the specific model defined by ^ = 1 + 



where 6 is a constant. This structure follows, for 
instance, by the assumption of invariancc with respect to 
a discrete Z2 symmetry under which the field ip is odd. 
The model has three free parameters: the initial scalar 
field configuration {ipjn,ip[^) and the constant b. 

It is well known that the system (13) has General Rel- 
ativity as a late time attractor solution [67]. Moreover, 
it has been found that, under very general assumptions, 
any departure from the GR attractor is associated to an 
enhancement of the Hubble rate [22, 68]. The strongest 
constraint on the model comes from BBN which states 
that the Hubble rate at the BBN can at most differ by 
10% from its GR predicted value [69-76]. We found, 
however, that the GR attractor is so efficient that even 
regions of the parameter space leading to pre-BBN val- 
ues of the ratio (15) order of magnitudes larger than one 
are allowed by the BBN bound. Fig. 20 shows the evo- 
lution of the ratio (15) for three benchmark points in 
the parameter space of our reference model: the solid 
line corresponds to the scalar field parameters fin = 1.9, 
= 0.45 and 5 = 8 (Model 1); the dashed line to 
ipin = 1.5, = 0.4 and = 8 (Model 2); the dot- 
dashed line stands for i^in = 1.5, (p[^ = 0.4 and 6 = 4 
(Model 3). 

The efficiency of the attractor is related to the evo- 
lution of the scalar field during radiation domination: 
throughout this epoch, when the temperature of the Uni- 
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FIG. 26: Enhancement factor 77* of the Hubble rate at the DM 
decoupling, for the three cosmological models of Fig. 20. The 
open circles, crosses and dots show the configurations which 
explain the PAMELA data without overproducing other sig- 
nals in terms of a thermal relic with the correct relic adun- 
dance and annihilation into e^e~ , fx^ ii~ and t'^t~ , respec- 
tively. The DM distribution is from Via Lactea II and the 
cosmic-ray propagation parameters are set at the MED case. 



verse crosses the value corresponding to the mass of a par- 
ticle in equilibrium in the thermal bath, the right hand 
side of the field equation becomes different from zero. 
This makes the field evolve towards its value at the GR 
attractor. Such a mechanism produces typical features in 
the field evolution (Fig. 21) and in the effective equation 
of state parameter (Fig. 22). 

For the scalar-tensor models introduced so far, we have 
then calculated the relic abundance for thermal relics 
able to explain the PAMELA data and compatible with 
the other astrophysical constraints, i.e. for the annihila- 
tion cross sections derived in Section II. Figs. 23, 24 and 
25 show the cases for pure leptonic annihilation channels: 
e^e~ , n'^ fJ-~ and r+r" respectively. The solid lines refer 
to Model 1, 2 and 3 introduced above and the dashed 
line shows the values of the relic abundance obtained in 
GR for the values of annihilation cross sections which 
are able to explain the PAMELA data within the astro- 
physical bounds. While in standard cosmology the relic 
abundance is typically too low (whence the need for a 
"boost"), for the three scalar-tensor comologies the relic 



FIG. 27: The same as in Fig. 26, for the cored isothermal 
halo profile. 

abundance is much larger and may be compatible with 
the WMAP value, represented by the horizontal line. We 
see that in the case of the specific Models 1, 2 and 3 pre- 
sented here, solutions to the PAMELA "anomaly" are 
found for specific values of the dark matter mass (in a 
range bewteen 100 and 500 GeV), dependent on the anni- 
hilation channel. These specific models therefore explic- 
itly implement a cosmological solution to the PAMELA 
"excess" in terms of a thermal dark matter particle with- 
out violating the observations discussed in Section II. 

Finally Fig. 26 and 27 show, in terms of the terminol- 
ogy of the previous Section on a generic deviation of the 
expansion rate from the GR behaviour, i.e. in terms of 
the enhancement parameter 77, the increase produced in 
Models 1, 2 and 3 here under discussion. Since, as it was 
discussed in Section IV, rj represents in our definition the 
value of the enhancement of the Hubble rate at a temper- 
ature corresponding to the freeze-out temperature T^q 
in GR, we have defined rj* shown in Fig. 26 as: 

l + V*^ AiT^^) = f^^(TpGR) (16) 

^GR 

We see that Models 1, 2 and 3 endow with enhancements 
of the order of 5 50. The open circles, crosses and 
filled circles in Fig. 26 (Via Lactea II) and Fig. 27 (cored 
isothermal) show the values of dark matter masses for 
which accomodate PAMELA positron data while being 
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compatible with the other bounds is possible, for the 
e+e~, and r+r" channels, respectively. We can 

see that for a cored isothermal halo profile, the cosmolog- 
ical model require slightely lighter dark matter in order 
to explain the PAMELA data, as compared to the Via 
Lactea II halo. 



VI. CONCLUSIONS 

Alternative cosmologies, based on extensions of Gen- 
eral Relativity, predict modified thermal histories in the 
Early Universe in the prc-BBN era, which is not directly 
constrained by cosmological observations. A typical pre- 
diction is that the expansion rate is enhanced with re- 
spect to the GR case: this, in turn, implies that ther- 
mal relics typically decouple earlier and with larger relic 
abundances. The correct value of the relic abundance 
for a thermal relic, value which cannot exceed the cos- 
mological determination of the dark matter content of 
the Universe, is therefore obtained for larger annihilation 
cross-sections, as compared to standard cosmology. 

Indirect detection rates of dark matter directly depend 
on the current value of {(Tannv) in the galactic halo. In 
the case of a dominant s-wave annihilation, which is typ- 
ical for most of the cold dark matter candidates in large 
portions of the parameter space of New Physics models, 
larger values of {a^nnv) required to match the WMAP 
dark matter abundance in modified cosmologies imply 
larger signals in the Galaxy. We have exploited this fea- 
ture in a twofold way. 

First of all, we can use the large host of independent 
results on the search for indirect signals of dark matter 
to set bounds on the enhancement of the Hubble rate in 
the pre-BBN era: this idea was introduced in Ref. [3], 
and there pursued by using exotic antiproton searches 
in cosmic-rays. A first attempt to use the gamma-ray 
signal was done in Ref. [4]. In the current Paper we ex- 
tend these analyses by introducing a whole set of indirect 
detection signals, which became increasingly relevant in 
the last months with the recent results from detectors 
like PAMELA, Fermi-LAT and HESS. We have classified 
cathegories of cosmological models with an enhanced ex- 
pansion rate and we have derived bounds on them under 
the hypothesis that the dark matter is a thermal relic. 
The observational data we have used to set bounds on 
the dark matter annihilation cross section (from which, 
in turn, we have derived the bounds on the cosmolog- 
ical models) come from radio observations, gamma-ray 
observations, inverse Compton photons, antiprotons and 
positrons in cosmic-rays and the optical depth of the Cos- 
mic Microwave Background photons. The bounds have 
been derived for a generic dark matter particle, by study- 



ing separately the different annihilation channels. 

The recent results on the measurement of the positron 
fraction realized by the PAMELA detector, have shown 
a clear and steady raise at energies above 10 GeV. This 
behaviour, which is currently under deep scrutiny, has 
one of its interpretations in terms of a dark matter signal 
through dark matter annihilation in our Galaxy. The 
theoretical analyses which have discussed this possibil- 
ity show that, in order to explain the PAMELA "ex- 
cess" , the annihilation cross sections need to be orders of 
magnitude larger than those required in standard cos- 
mology to explain the observed amount of cold dark 
matter. This fact poses a problem for a thermal relic, 
and various mechanisms have been invoked to boost the 
positron signal without spoiling the correct value of relic 
abundance (astrophysical boosts, Sommerfeld enhanced 
cross-sections). Since in cosmologies with an enhanced 
expansion rate we naturally require larger annihilation 
cross-sections, we have discussed the properties of these 
alternative cosmologies in order to be able to explain the 
PAMELA "puzzle" without violating other observations. 
We have derived the required amount of enhancement for 
different cosmologies. 

Finally, in order to demonstrate the feasibility of such 
an approach, we have discussed scalar-tensor theories of 
gravity, for which we constructed explicit models able to 
reproduce the required cosmological boost to explain the 
PAMELA data while being compatible with astrophysi- 
cal constraints. 
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